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ABSTRACT 


In-flight measurements of the SR-71 near- field sonic boom were obtained by an F-16XL airplane at flightpath 
separation distances from 40 to 740 ft. Twenty-two signatures were obtained from Mach 1 .60 to Mach 1 .84 and 
altitudes from 47,600 to 49, 150 ft. The shock wave signatures were measured by the total and static sensors on the 
F-16XL noseboom. These near-field signature measurements were distorted by pneumatic attenuation in the pitot- 
static system. The near-field pressure signatures were corrected by modeling the magnitude and phase of the pitot- 
static sensors and accounting for their effects using optimal deconvolution. Measurement system magnitude and 
phase characteristics were determined from ground-based step-response tests and extrapolated to flight conditions 
using analytical models. Deconvolution was implemented using Fourier transform methods. Comparisons of the 
shock wave signatures reconstructed from the total and static pressure data are presented. The good agreement 
achieved gives confidence of the quality of the reconstruction analysis. Although originally developed to reconstruct 
the sonic boom signatures from SR-71 sonic boom flight tests, the methods presented here generally apply to other 
types of highly attenuated or distorted pneumatic measurements. 

INTRODUCTION 

When an aircraft flies faster than the local sonic velocity, pressure discontinuities develop at the vehicle surface 
and propagate outward to form shock waves. These waves form in places where sharp changes in the vehicle 
geometry occur, such as at the nose, tail, canopy, wing and body junctions, and engine inlets. The flow-field region 
where shock waves exist for each geometric discontinuity extends an approximately 250-ft radius from the aircraft 
and is referred to here as the near-field region. With increasing distance from the aircraft, variations in the local 
speed of sound cause the various shocks induced by geometric changes to gradually coalesce into two major shock 
waves (refs. I, 2). This coalescence gives rise to the characteristic N-wave profile which sweeps along the ground 
and produces the well-known sonic boom effect. The region where the shock waves have coalesced into a single 
bow shock and a single tail shock is referred to here as th tfar-field region. The region where coalescence has begun 
but is incomplete is referred to as the mid-field region. Figure 1 shows the three sonic boom regions for the SR-71 
aircraft (Lockheed Aerospace, Burbank, California). 



Figure 1 . Near-, mid-, and far-field sonic boom signatures. 
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The perceived strength of a sonic boom is proportional to the peak overpressures which occur in the far-field 
region. Atmospheric factors which cause booms of significant strength to dissipate are not well understood. 
Analytical computations which describe the dissipative forces in the atmosphere were performed, but these 
computations are not well-supported by correlation with a broad base of experimental flight data (ref. 3). To 
benchmark the analytical code predictions, correlation of ground-based, far-field, sonic boom signature 
measurements with those measured in the near and midfields is necessary. 

Data are relatively easy to obtain in the far-field, and a significant database exists for far-field measurements. To 
complete the picture, quality near-field sonic boom signatures are required. Unfortunately, because of the 
complexity of obtaining controlled measurements while operating in the near field of a supersonic aircraft, research 
quality signature measurements in the near field have never been satisfactorily obtained in-flight. To demonstrate 
the feasibility of making near-field measurements along the vertical centerline of a large supersonic aircraft at 
supersonic speeds, an experiment was conducted at NASA Dryden Flight Research Center (DFRC), Edwards, 
California (ref. 4). A probe aircraft was used to sense the near-field shock wave signature of an SR-71 aircraft. For 
these tests, an F-16XL aircraft (General Dynamics, Fort Worth, Texas) was used as the probe aircraft. The F-16XL 
aircraft was chosen because of its high Mach number flight envelope and its supersonic endurance. Shock wave 
signatures were measured by the F-16XL noseboom pitot-static system. 

Conceptually, pitot-static pressure time histories can be correlated with vehicle trajectory data to infer the spatial 
distribution of the near-field signature. In practice, however, the task is not as straight forward. Raw measurements 
resulting from the F- 1 6XL pitot-static system were badly distorted. These pressure distortions were primarily caused 
by pneumatic attenuation in the airdata measurement system. In addition, the local flow-field of the F-16XL aircraft 
affected the raw pressure data. 

This paper describes a novel method for reconstructing the sonic boom signature by modeling the magnitude and 
phase characteristics of the pitot-static systems and accounting for their effects using optimal deconvolution. 
Pressure sensor magnitude and phase were determined from ground-based step-response tests and extrapolated to 
flight conditions using analytical models. To remove the F-16XL local flow-field effects, reference total and static 
pressure time histories were computed. These reference time histories represent the total and static pressure fields 
which would occur in the absence of the SR-71 aircraft. The reference time histories were constructed from data 
sources independent of the F-16XL pitot-static measurement system. To evaluate the shock wave signatures, the 
reference pressure time histories were subtracted from the deconvolved pressure time histories. Examples of 
reconstructed flight signatures are presented in the Flight Test Procedures, Analysis, and Results section. 

BACKGROUND 


Mullens (ref, 5), Smith (ref, 6), and Maglieri (ref. 7) have obtained in-flight signature measurements. Mullens’ 
data were obtained from 85 to 1770 ft below and from 100 to 1420 ft to the side of an F-100 airplane at Mach 1 .05. 

A General Dynamics (Fort Worth, Texas) B-58 bomber airplane served as the boom generator during the Smith 
and Maglieri investigations. Smith’s results include data for the F-100 and F-104 airplanes at Mach 1 .2 and data for 
the B-58 from Mach 1.3 to Mach 1.8. Data were obtained for lateral separation distances from 120 to 425 ft. The 
investigation was conducted to determine the near-field-flow patterns for assessing the dangers of close- formation 
flying. 

Maglieri’s data were obtained at distances from 1300 to 9100 ft below and from 1700 to 2000 ft above the B-58 
aircraft from Mach 1 .42 to Mach 1 .69. Emphasis was on obtaining data to assess the lift and volume interaction 
effects on the near-field signatures of a delta- wing bomber. For the experiments conducted in refs. 5-7, the data arc 
either too low in Mach number, not below the generating aircraft, or too far away for the present needs of bench- 
marking analytical computations. These results are summarized in table 1. 
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Table 1 . Previous sonic boom signature measurement flight tests (refs* 5-7). 


Investigator 

Aircraft 

Mach 

number 

Distance 
and position 

Mullens, 1 956 

F-100 

1.05 

85 -1770 ft below 

100 -1420 ft to side 

Smith, 1960 

B-58, F-100, F-104 

to 

1 

oo 

120-425 ft to side 


Maglieri, 1963 

B-58 

1.42-1.69 

1300 - 9 100 ft below 

1 700 - 2000 ft above 


VEHICLE AND INSTRUMENTATION DESCRIPTIONS 

As mentioned in the Introduction section, the SR-7 1 airplane was chosen as the shock-generating vehicle, and the 
F-16XL airplane was chosen as the probe aircraft. This section describes these vehicles and discusses the 
instrumentation systems available on each vehicle. 

SR-71 Aircraft 

Figure 2 shows the SR-71 aircraft. This aircraft is one of three at DFRC that is on loan from the U. S. Air Force. This 
two-place, delta-wing aircraft was formerly used for military reconnaissance. Powered by two J-58 axial-flow turbojets 
(Pratt & Whitney, West Palm Beach, Florida), an SR-71 aircraft can sustain a cruise of Mach 3.2 at an altitude of 
85,000 ft (ref. 8). The vehicle instrumentation used in this analysis includes the inertial navigation system (INS), astro- 
navigation system (ANS), and airdata measurements. These airdata measurements are teed off from the production 
airdata probe to NASA-supplied transducers. Data from these sensors were intercepted from the ship communication 
bus by a pulse code modulation (PCM) system interface and telemetered to ground for real-time monitoring and 
postflight analysis. Flight trajectory was monitored and recorded by ground-based, C-band radar. The radar data stream 
was merged with the aircraft telemetry data stream postflight for research analysis. 



EC 92 09241-1 


Figure 2. The SR-71 aircraft. 
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F-16XL Aircraft 


Figure 3 shows the F-16XL airplane. This airplane is powered by a nonproduction F11Q-GE-100 turbofan engine 
(General Electric, Lynn, Massachusetts) (ref. 9). The F-16XL aircraft has been flight tested to a maximum of 
Mach 2.05 and a peak altitude of 65,000 ft. The vehicle is outfitted with a research-quality instrumentation system with 
airdata being obtained with a nose-boom-mounted pitot-static airdata probe. Total pressure measurements were 
provided by a single pitot port at the boom tip, and static measurements were provided by a partial ring of static pressure 
ports located circumferentially on the boom. For this analysis, total and static pressure were teed off from the noseboom 
pressure supply lines and were sensed by a set of identical absolute digital pressure transducers whose full-scale ranges 
were 0-2736 psf (0-19 psi) with a digital resolution of 18 bits (0.01 psf). The total and static pressure transducers 
had measured 3 standard deviation, a, accuracies of approximately ±0.0075 percent of the calibrated full-scale 
(± 0.21 psf). In addition to the research instrumentation, the total and static pressure inputs supplied the ship’s system 
central airdata (CAD) computer and the cockpit instrumentation gauges. This additional instrumentation added 
considerably to the volume entrapped within the sensing systems and contributed largely to the measurement 
magnitude and phase distortions described in the Ground Test Procedures, Analysis, and Results section. The total 
volume entrapped in the static-pressure-measuring system was roughly four times that of the total pressure system. 

Other research measurements included inertial navigation parameters from the INS and space positioning and 
velocity data from a ground-based radar track of the C-band beacon. The radar data stream was merged with the aircraft 
telemetry data stream postflight for research analysis. Onboard data measurements were collected by the 10-bit PCM 
system, recorded onboard with an analog tape recorder, and telemetered to ground in encoded form. The 1 8-bit, airdata 
pressure measurements were collected as two 10-bit words. These words were concatenated on the ground during 
postflight data analysis. 
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RESULTS AND DISCUSSION 


This section summarizes the results of ground-based, step-response tests used to determine the pressure sensor 
magnitude and phase characteristics, flight test maneuvers used to obtain the raw pressure data, and analysis used to 
reconstruct the sonic boom signatures. Results from the flight tests, including representative boom signature data, are 
presented. 

Ground Test Procedures, Analysis, and Results 

This subsection describes the ground test procedures used to estimate the airdata system parameters. Analysis 
methods are described, and results of the ground tests are presented. 

Procedures 

Because spatial distribution of the near-field, sonic boom signature must be inferred from the time correlation of 
pressure and space positioning data, the magnitude and phase characteristics inherent to the F-16XL pitot-static system 
must be well-known. Unfortunately, the internal pneumatic layout of the CAD computer and the cockpit displays is 
unknown, and accurately modeling the responses of the system components analytically is impossible. To determine 
the baseline response characteristics of the measurement system, a series of simple, ground-based, step-response tests 
was performed. 

The step-response test procedures are described next. First, a pressure fitting was attached to the total pressure 
portion of the noseboom. Then, a “Y” connector was attached to this fitting. One branch of the “Y” was connected to 
an accurate reference pressure calibration unit. The other branch of the “Y” had a short length of tubing called the vent 
line. The vent line was doubled over on itself and held sealed with a surgical hemostat. Next, the total pressure side of 
the airdata system was then evacuated or pressurized to the desired level above or below the ambient pressure. When 
the pressure stabilized, the tubing between the calibration unit and the noseboom was doubled over and sealed with a 
second hemostat as close as possible to the noseboom. To vent the system to ambient pressure condition, the first 
hemostat on the vent line was rapidly opened allowing a sharp pressure change to be introduced to the total pressure 
system. After the system reached ambient pressure level, the vent line was resealed using the hemostat and pressurized 
(or evacuated) to the next test condition. Lastly, after the total pressure tests were completed, step-response tests on the 
static pressure system were conducted using an identical procedure. 

To evaluate the range of pressure rises, ground tests were performed for positive and negative step inputs with data 
being acquired for steps with ±0.1, ±0.5, and ±1.0 psi magnitudes. These data responses were recorded using the 
F-16XL onboard data acquisition system at a rate of 25 samples/sec. Ambient conditions were recorded off-line using 
a hand-held barometer and thermometer. These data were used later in the posttest analysis. For each condition, the 
step-response tests were repeated twice. 

Analysis 

For small step inputs, laminar flow conditions exist in the pneumatic system. The analysis presented in appendix A 
demonstrates that total and static airdata system responses can be closely approximated by a second-order system. 
The transducer measurements are the system output, and external pressure provides the system input. Appendix A 
also presents a method to extrapolate the ground test natural frequencies and damping ratios to flight conditions. 
Appendix B describes the parameter identification method used to extract the system parameters from the ground tests. 

Results 

The results of the step-response tests were averaged for the total and static pressure systems. Table 2 lists the 
averaged natural frequencies, damping ratios, and time lags. 

Selected comparisons of the second-order response fits to the measured data are presented in figures 4(a) and 4(b) 
for the total pressure system and in figures 5(a) and 5(b) for the static pressure system. Clearly, the second-order model 
provides a good fit to the ground-based data. 
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Table 2* Ground test results for the second-order system parameters. 


Pressure 

Natural radian 
frequency, 

Damping 

Time 

system 

rad/sec (Hz) 

ratio 

lag, sec 

Total 

17.47 (2,780) 

0.841 

0.096 

Static 

8.485(1.351) 

2.254 

0.531 


Normalized 

pressure 



Normalized 

pressure 



(b) Sp = — 0.493 psi. 


Figure 4. Normalized step-response tests: Comparisons of model fits with data for total pressure system. 
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Normalized 

pressure 



Normalized 

pressure 



(b) 5 p = -0.495 psi. 

Figure 5. Normalized step-response tests: Comparisons of model fits with data for static pressure. 


Flight Test Procedures, Analysis, and Results 

This subsection describes the flight test procedures used to gather the raw sonic boom data. Methods used to 
reconstruct the sonic boom signatures are described, and results of the analysis are presented. Comparisons of boom 
signatures reconstructed using static and total pressure data are presented. 

Procedures 

The F-16XL aircraft was positioned on the vertical plane of symmetry behind and below the SR-71 aircraft at the 
required vertical separation distance. As the F- 1 6XL aircraft probed forward into the shock field of the SR-7 1 aircraft, 
the pilot could see when the F- 16XL passed through the shock waves. The F-16XL pilot described this shock as looking 
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“like a sheet of water.” The encounter times were marked for later correlation with telemetered data. In all, 
22 signatures were obtained from Mach 1.60 to Mach 1.84 and at altitudes from 47,600 to 49,150 ft. Flightpath 
separation distances between the vehicles varied from 40 to 740 ft. 

The desired overpressure signature is the difference between the corrected pressures sensed at the noseboom and 
a reference pressure. Reference pressure represents the pressure that would have occurred in the absence of the 
SR -71 aircraft. As the F-16XL aircraft probed the SR-71 shock field, interactions with SR-71 shock waves made it 
impossible for the probe aircraft to maintain a steady separation and velocity through the flow field. These speed 
and altitude variations significantly affected the pressure values sensed at the noseboom. Thus, the reference 
pressure data must be reconstructed using trajectory data that are measured independently from the F-16XL 
noseboom. This reconstruction is discussed in the Analysis subsection. 

This procedure was used to reconstruct the respective sonic boom signatures. Results of the ground-based, step- 
response tests were extrapolated to flight conditions using the method described in appendix A. The raw pitot-static 
pressure measurements were corrected for the effects of pneumatic distortion using the method described in 
appendix C. Next, the reference pressure time histories were generated and subtracted from the deconvolved 
signatures to give the shock wave overpressure measurement. For this analysis, signatures from the static and total 
pressure measurements were analyzed separately. 

Analysis 

Data sources used to generate the reference pressure data were radar space positioning and velocities (ref. 10), INS 
attitudes and velocities, and weather parameters. These weather parameters included windspeed and direction, 
ambient pressure, and ambient temperature from meteorological analysis or rawindsonde weather balloon data and 
National Weather Service charts (ref. 11). The static pressure measured by the static pressure sensor before 
penetration of the SR-71 shock wave field was used as the baseline zero overpressure. As the F-16XL aircraft probed 
the shock field of the SR-71 aircraft, the reference static pressure data were generated using variations in radar- 
derived altitude and the hydrostatic equation (ref. 12). To give variations in the total pressure baseline, a reference 
airspeed was computed from inertially derived ground speed data (either radar or INS), and windspeeds were 
computed from the meteorological analysis. Temperature from the meteorological analysis was used to compute the 
reference Mach number from the reference airspeed data. Using the Rayleigh pitot equation (ref. 14), the reference 
total pressure data were generated using the reference Mach number and static pressure data. 

Results 

Sample sonic boom time history signatures are presented next. These data are representative of all results obtained 
during the sonic boom flight tests. Comparisons of the raw and deconvolved pressure signatures is presented first. 
Next, the reference pressure data are subtracted to yield sonic boom overpressure time histories. Lastly, signatures 
obtained using total and static pressure data are evaluated independently, and the results are compared. 

Figure 6 illustrates the effect of the deconvolution on the raw signature measurements. These data were obtained 
in the near field at a vertical separation distance of 200 to 260 ft below the centerline of the SR-71 aircraft. Here, the 
shock wave field was probed from the rear of the SR-71 aircraft moving forward, Time histories of raw and 
deconvolved static and total pressure are shown. Shock waves generated by the aircraft tail, inlet, canopy, and bow 
are labeled. Raw total pressure is slightly lagged from the deconvolved data and has a slightly higher amplitude. In 
this case, the total pressure sensor is mildly resonant, and the deconvolution algorithm acts as a signal attenuator. 
On the other hand, when compared to the deconvolved static pressure data, the raw static pressure has been greatly 
lagged and attenuated. The raw static pressure is more heavily attenuated than the total pressure measurement 
because the system pressure level is considerably lower. In addition, the entrapped volume in the static pressure 
system is larger by a factor of approximately four. See appendix A. 
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Raw 




Figure 6, Comparisons of raw and deconvolved shock wave pressure signatures (200 - 260 ft below vertical separation). 


The reference pressure data represent the pressures that would be sensed in the absence of the SR-71 aircraft. Using 
the data of figure 6, comparisons of the deconvolved pressure data and the reference data for the static and total 
pressures are presented in figures 7(a) and 7(b), Note that little variation exists in the reference static pressure. Because 
of probe aircraft speed and separation variations as the F-16XL aircraft penetrated the SR-71 shock field, large 
variations occurred in reference total pressure during the maneuver. 

Using the data of figures 6 and 7, comparisons of signatures derived from the deconvolved static and total pressure 
data are presented in figure 8. Raw signatures reconstructed by subtracting the reference pressures from the 
undeconvolved pressure measurements are also presented. Clearly, the raw signatures are radically different for the 
total and static pressures. On the other hand, the deconvolved total and static pressure signatures show close agreement. 

Note that the static and total pressure signatures were obtained from fundamentally different How-field 
measurements. In addition, at high Mach numbers, shock wave losses cause the overpressures measured by the two 
types of sensors to disagree. However, for these test conditions, simple adiabatic analysis (ref. 13) showed that in the 
absence of other measurement errors overpressures measured by the static and total pressure measurements should 
closely agree. Because these total and static pressure signatures were derived independently, the good agreement 
between the signatures gives confidence in the validity of the deconvolution and trajectory analysis. Although 
deconvolution analysis gives a good representation of overpressure time histories, because of jostling of the probe 
aircraft as it penetrated the SR-71 shock field, these data were not obtained at a constant separation distance. 

Figure 9 shows similar comparisons for the midfield (580 -740 ft below). For these data, the field was measured from 
front to aft as the probe aircraft decelerated. Thus, the signature appears to be backward when compared to figure 8. 
Notice in figure 9 that there is no distinct canopy shock. At this large separation distance (580-740 ft), the canopy 
shock has been absorbed by the inlet shock. This absorption is a normal feature of the shock coalescence process. 

Finally, signatures measured in the very near field (40-90 ft below) are presented in figure 10. For this maneuver, 
the field probe was performed below and to the side of the SR-71 aircraft. This causes the engine inlet shocks to appear 
as two distinct shocks. For the near-field data of figure 10, the very large tail shock strength is underpredicted by the 
static pressure data. On the other hand, for the signatures measured in the mid- and far-fields, the static and total 
pressure data give approximately the same tail shock strengths. Analysis presented in appendix D suggests that this 
discrepancy results from turbulent flow conditions being induced in the static pressure system by large pressure 
gradients. When turbulent flow is induced in the static pressure system, the tube damping rises dramatically. In addition, 
the second-order deconvolution model, which uses laminar flow assumptions, underpredicts the size of the shock wave. 
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(a) Static pressure. 



(b) Total pressure. 

Figure 7. Comparisons of deconvolved and reference pressure data (200 - 260 ft below vertical separation distance). 
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Figure 8. Comparisons of static and total pressure derived sonic boom signatures (200 - 260 ft below vertical 
separation distance). 



Figure 9. Midfield shock overpressure signature comparisons (580 - 740 ft below vertical separation distance). Note 
that the field was probed from front to back. 
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Figure 10. Comparisons of sonic boom signatures for very near field (40 - 90 ft below vertical separation distance). 

Figure 1 1 presents the approximate boundaries for the maximum allowable difference between the external and 
measured pressures (Pq and P^) which can be tolerated without inducing transitional flow within the tube. Here, the 
shaded regions will have maximum Reynolds numbers which are less than the tube flow transition Reynolds number of 
2300. Total pressure can have higher step inputs with increasing Mach number without inducing transitional flow. At 
the flight conditions ol figures 6 through 9, the total pressure system (fig, 11(a)) can sustain a maximum pressure 
difference in the tube of nearly 50 psf without going turbulent. The static pressure system (fig. 1 1 (b)) can sustain a 
maximum difference of approximately 40 to 45 psf without becoming turbulent. 

For the near-field signature data, turbulent flow is induced in the static pressure system by the tail shock because the 
system response latency allows a large pressure gradient to form. This pressure gradient induces sizable tube flow 
velocities, and the critical transition Reynolds number is exceeded (ref. 18). Turbulent flow does not form in the total 
pressure-measuring system because the system response equalizes the pressure within the tube before large flow 
velocities can build up. For the very near-field data of figure 10, the differences between the deconvolved total pressure 
signature and raw total and static pressure signatures are plotted in figure 1 2. Clearly, the pressure differential level for 
the total pressure system is well below the level required to reach critical Reynolds number. On the other hand, at the 
tail shock, the static pressure differential exceeds the critical threshold. 




(a) Total pressure system. (b) Static pressure system. 

Figure 11. Maximum allowable pressure differential in systems for laminar flow. 
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Figure 1 2. Maximum internal pressure differential for very near-field total and static pressure measurement systems 
(40 to 90 ft below vertical separation distance). 


CONCLUDING REMARKS 

In-flight measurements of the near-field SR-71 shock wave pressure field were obtained at separation distances 
ranging from 40 to 740 ft by probing with the noseboom pitot-static system of an F-16XL airplane. The raw pitot- 
static measurements were distorted by pneumatic distortion and influences of the probe aircraft flow field. This 
paper develops a novel method to correct for these distortions. 

Pneumatic distortions were corrected by modeling the pneumatically induced magnitude and phase distortions and 
by accounting for their effects using optimal deconvolution. The local flow-field effects were removed by 
construction reference time histories of the total and static pressure which would occur in the absence of the SR-71 
shock waves and subtracting the results from the deconvolved pressure data. The pitot-static measurement system 
magnitude and phase characteristics were determined from ground-based step-response tests and extrapolated to 
flight conditions using analytical models. The deconvolution is implemented using Fourier transform methods. 
Reference time histories for total and static pressure were constructed using data sources independent from the pitot- 
static measurement system, including C-band radar tracking, inertial navigation system attitudes and velocities, and 
rawindsonde weather balloon sounding data. 

Twenty-two signatures were obtained from Mach 1.60 to Mach 1.84 and altitudes from 47,600 to 49,150 ft. 
Examples of reconstructed flight signatures were presented. Agreements of signature shape and overpressure 
magnitude for the static and total pressure-derived signatures are excellent. This agreement gives confidence of the 
quality of the reconstruction analysis. 

Although developed to reconstruct the sonic boom signatures from the SR-71 sonic boom flight tests, the methods 
described here apply to other types of highly attenuated or distorted pneumatic measurements. Other applications 
might include the reconstruction of pressure transients in piping systems from single point measurements and 
fidelity enhancement of unsteady pressure measurements. 


Dryden Flight Research Center 

National Aeronautics and Space Administration 

Edwards, California, April 5, 1995 
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APPENDIX A 


PITOT-STATIC SYSTEM RESPONSE MODEL 

This appendix presents the pneumatic system model and reduces the model to a quasi-second-order system. Rules 
for extrapolating the system parameters to various flight conditions are developed for laminar flow conditions. 

System Model 

For simple tubing geometry (fig. A-l), the frequency response characteristics of the system may be expressed in 
closed form as (ref. 14) 


w 


cosh JaL/c j + sinh Jo-L/ c 


where a is the square of the damped radian frequency 


2 R a 


.2. - R a 


a = s + — s = - to + j — (a=> s= id) 

P P 


where 


P 

P 


L 

0 


S = jo 


L 


c 

D 


measured pressure 
input pressure 
radian frequency 
tube length 
local speed of sound 
tube diameter 


(A-l) 


<A-2) 


where R a is the acoustical impedance of the system, and A c 


kD /4 is the cross-sectional area of the tubin 


£T 


P L 


\ 

...... ’ 

D J — u(r,x,t) 

.. , , .. , . . . ; ; 

, .. ; : ' .. " ; ; ; ; ; ; ; 

t 

j.» - ' 

' J..'. ■ ' 1 ■ '.O ■ "> «. txgf 



, . . .. 

L ► 

950423 


Figure A-l . Simple tubing geometry. 
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The acoustical impedance of the system is the resistance to flow with the tube and is proportional to the surface 
skin friction coefficient and the Reynolds number (refs. 15, 16), 


where 


*«- 


2 V R ey C f 

D 2 


c 

f = skin friction coefficient 

ft 

ey = Reynolds number 

It = dynamic viscosity 


(A-3) 


For one-dimensional tube flow with fully developed laminar flow conditions. 


c 

C f~R 

ey 


(A-4) 


Thus, 


if =32p 
* D 2 


(A-5) 


System Model Reduction 

For complex plumbing geometries in which branches or cascades occur, the basic solution, equation A-l, may be 
used as a building block to develop end-to-end solutions by modifying the boundary condition to allow recursive 
implementation. The frequency response at each node is given as a function of the frequency response of the 
previous node. This method, however, requires accurate knowledge of the system configuration. Unfortunately, the 
F-16XL airdata system has a complex plumbing arrangement with multiple lines, tees, compression fittings, trapped 
volumes, and transducer connections. As a result, knowing the precise geometric configuration of the system is 
difficult. Thus, analytical computations based on geometry assumptions would be highly suspect. Instead, it is 
assumed that for the sonic boom data set, the F-16XL systems can be represented by equivalent “single node” 
systems in which the damping and frequency are determined from ground experiments and extrapolated to flight 
conditions. Here, the geometry effects are determined using the ground-based, step-response data, and the results 
are extrapolated to flight conditions using physical rules. 

To perform the analysis, first, the form of the analytical model must be reduced to a manageable form, so 
extrapolation rules can be developed. This model reduction is done by performing a Taylor series expansion of 
equation A- 1 , where 


cosh 


JaL 

c 



m - 0 



(2m)! 


(A-6) 
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and 


sinh 


JaL 


L~\ 2m+l a (2m+l)/2 


- E 

m = 0 


[;] 


(2m + 1)! 


Substituting into equation A-1 and collecting like powers of a 
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Evaluating equation A-2 where p is the longitudinal average of the tubing density and regrouping terms 
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If all terms (in s) of equation A-9 of order greater than 2 are neglected, then a reduced order model which describes 
the behavior dominant harmonic results, 
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For fully developed laminar flow conditions, the acoustical impedance is a constant. In addition, equation A- 10 is 
a linear filter of the form 
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with natural frequency, damping ratio, and time lag 
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Development of Extrapolation Rules for Fully Developed Laminar Flow Conditions 

Clearly, equation A-14 contains two groups of terms: those which depend on the geometry only and those which 
depend on properties of the local flow, such as p, p, and c . Because the geometry is fixed for a given measurement 
configuration, the time constant may be extrapolated to various flight conditions (assuming constant y) according 
to the rule 
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Here again, the values for p , p , T, and p are the longitudinal averages. The expression for natural frequency 
(equation (A- 12)) does not neatly decouple into flow and geometry terms; 
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For laminar flow conditions, however, it can be shown with an order of magnitude analysis that 
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and the extrapolation rule reduces to the practical relationship 
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The damping ratio is computed directly from the natural frequency and time lag. The second-order filter given by 
equation A- 11 with parameters given by equations A- 12 and A- 13 will be used to perform the deconvolution 
analysis. 
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APPENDIX B 


ESTIMATION OF SYSTEM PARAMETERS FROM STEP-RESPONSE DATA 


The system parameters, % and (a n , were extracted from the step-response data using a linearized least-squares 
technique. As discussed in appendix A, the system response is described as a second-order system whose frequency 
domain equation can be written as 


P L (s) = 
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(B-l) 


Where P L is the measured pressure, P Q is the input pressure. 

Transformation to Discrete Filter 

If the second-order system model is discretized using the bi-linear transform (ref. 16), 
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Where, AT is the sample interval of the system. Then when equation B-2 is substituted into equation B-l and like 
powers in z are collected, the resulting discrete transfer function is 
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Transforming to the time domain using the transformation rale yields (ret. 16) 


z"[*(z)l = x k + n 


(B-8) 


The resulting difference equation is 

mv Lt ^ A mf Li ^m\p Lk = (B ' 9) 

with k as the time index. 

Formulation of Matrix Equations and Least-Squares Solution 

If TV + 1 data points are available, then equation B-9 may be written as the over-determined system of linear 
equations. For example. 
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Using the pseudoinverse method (ref. 1 8), which is equivalent to performing linear least-squares estimation* to 
solve equation B-10 for the difference equation coefficients, A[ 0] * A[ 1 ] , A[2\ , then the natural frequency and 
damping ratio can be computed by manipulation of the original relationships 
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APPENDIX C 


DEVELOPMENT OF OPTIMAL DECONVOLUTION ALGORITHM 

In this appendix, the algorithm for reconstructing the “true” input pressure from the raw data is developed. This 
algorithm is based on the methods of spectral convolution, and the resulting algorithm is implemented as a frequency 
domain filter. Conceptually, it is a trivial matter to compensate for the effects of the acoustical distortion by taking 
the numerical Fourier transform of the raw pressure data, multiplying by the inverse of the transfer function, and 
taking the numerical inverse transform. This method is referred to as spectral deconvolution. In practice, however, 
the implementation is not easy. 

Clearly, the measurement system will have some nominal signal-to-noise ratio, Rs n , which represents the 
fundamental accuracy of the measurement system. The noise is introduced into the measured signal by the 
transducer, its associated electronics, and the quantization of the measurement system. The measurement noise is 
typically white over the range of the sampled spectrum; however, this noise may contain concentrated regions oi 
spectral energy, such as noise introduced by electromagnetic interference. 

Regardless of the noise characteristics of the system, for pneumatic data at higher frequencies, the measured signal 
is attenuated. In addition, the signal-to-noise ratio drops considerably below the nominal value. At high frequencies 
where the measured signal is attenuated, the deconvolution transfer function acts as an amplifier. This function will 
amplify the attenuated signal and the additive measurement noise. This overamplification produces a very noisy 
reconstructed signal. Thus, an inversion filter which controls and minimizes noise amplification must be developed. 

Development of Compensation Filter 

Assuming that the dynamics of the pneumatic distortion are described by the second-order model developed in 


where pj is the measured pressure value, and p 0 is the input surface pressure. Then, the noisy pressure measure- 
ments can be modeled by the convolution 


where u(r - x) is the convolution function, and q(r) is additive measurement noise. For this analysis, it is assumed 
that the measurement noise is unbiased. That is, JSfq (t)\ = 0 . Because the noise is introduced by the measurement 
system, it is uncorrelated with the surface pressure. For example, E[p 0 (OTl(*)l = 0 , 

The Fourier transform of equation C-2 is 
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where P L ( to), P 0 (©), and //(to) are the Fourier transforms of the measurement, input, and noise function time se- 
ries, and © is the radian frequency parameter. Because the noise in unbiased, 


E[H(c 0)1 = £{F[rKr)]} = [ e m E[j\{t)]dt = 0 

* —OQ 


(C-4) 


Similarly, because the input pressure and measurement noise are uncorrelated, 
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The goal is to select an inversion filter, G(©), such that the surface pressure estimate minimizes the squared 
difference between the estimated surface pressure and the actual value; that is, seek an optimal filter of the form 


P 0 {©) = G(©)P L (©) 


(C-6) 


which minimizes the cost functional 
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subject to the constraints of equation C-l . The optimization method used to develop this filtering algorithm was first 
developed by Norbert Weiner. The resulting filter can correctly be referred to as a Weiner Filter (ref. 19). 
Substituting for Pq in equation C-7 
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Dropping the (©) notation and expanding terms in equation C-8, the cost functional becomes 
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This notation is assumed implicit for simplicity of notation. Defining 
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Letting Real(G) = G r , lm(G) = G im , Real(A L0 ) = A LOf , Im(A i0 ) - A L0 ^ ; then, G G r + jG im , 

A r s A , n + ; A T n .In addition, and substitution into equation C- 1 0 gives 
LO LU r LU i m 
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= (G r 2 + Gj)A iL -2(G/ i0r + G j . m A i0 J + A 00 (C-U) 

Because values for A IO , A LI , and A ()0 are defined by the expectation values for the measured and true input 
pressures (defined above), G r and G im are the only free parameters of the system, and necessary and sufficient 
conditions for optimization (ref. 1 5) of the cost index are (holding frequency constant) 
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Computing the partial derivatives yields 
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The optimal estimate of the input pressure spectrum based on the measured pressure spectrum and the constraints 
of the dynamics model is 
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The numerator of equation C- 14 depends on the correlation of the measured pressure to the true input pressure, an 
unknown value. Using the measurement model, equation C-6, however, equation C-14 can be written in terms of the 
expected values for the input pressure and the measurement noise. 
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Defining ||P 0 || 2 = E[P Q *(G>)P 0 (a>)] , Iffll* = the optimal deconvolution algorithm becomes 
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Equation C-17 weights the Fourier coefficients of the deconvolution as a function of the measured signal-to-noise 
power spectrum ratio, 

l f 0< m >ll (C-18) 

l|H(co )|| 2 

which is closely approximated by the nominal signal-to-noise ratio of the system. Clearly, as the signal-to-noise ratio 
increases, the optimal deconvolution algorithm becomes a simple deconvolution. 
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Equation C-17 is the final form of the optimal deconvolution algorithm. 

Evaluation of Signal-to-Noise Ratio for the SR-71 Sonic Boom Tests 

For the sonic boom analysis, the measurement noise was assumed to be white. Based on the manufacturer s 
accuracy specifications of ±0.0075 percent of full-scale, and a full-scale value of 2736 psf (19 psi), the nominal 
signal-to-noise ratio of the measurement system is 
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Here, , is the mean square of the measured time history over the course of the maneuver. 


Effects of Time Series Truncation 


For practical implementation using the fast Fourier transform (FFT), great care must be taken to avoid finite-sized 
windowing effects (ref. 19). Such effects would be amplified by the deconvolution algorithm. Truncation of a time 
series at t ± t c is equivalent to multiplying by a “truncation boxcar.” The “resolution bandwidth, the smallest 
frequency which can be independently resolved by the finite length Fourier transform, is inversely proportional to 
the length of the time window. When power is distributed at frequencies which are separated by less than the 
resolution, these bandwidths are smeared together. For this analysis the effects of truncation smearing can be 
mitigated by using a “taper window” before the transform of the original time series. Window selection has been 
researched extensively (ref. 19). Generally, the focus has been on choosing a truncation function so that the effects 
induced by the frequency domain convolution are reduced. For this analysis, a “cosine taper” window is used to 
mitigate these effects* This window is defined as 





(C-21) 


Here, k is a parameter which gives a degree of freedom to the window, and it is selected by engineering judgment. 
For the sonic boom analysis, a value of k = 5 was sufficient* 

Implementation of the Inverse Transform 


To ensure that a real-time series results when the frequency spectrum is inversely transformed, the discrete 
spectrum should be skewed symmetric about the Nyquist frequency* For the discrete series, the upper half of the 
spectrum should obey the relationship 
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As with the forward Fourier transform, the inverse transform also has windowing effects. The inverse transform 
is approximated by a finite sum. The truncated harmonics, which are a “window in the frequency domain, induce 
a convolution in the time domain. This convolution causes unwanted oscillations in the resulting time series, a 
process known as Gibb's phenomenon (ret . 19). To eliminate these oscillations, the deconvolved time series is 
filtered by taking weighted averages of adjacent time points, such as 
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The parameter, X , is a degree of freedom that must be selected using engineering judgment. The filter weights are 
arbitrary as long as they satisfy the relationship 


to prevent a change in the direct current (dc) level of the signal. For this analysis, a window length of 9 { X - 4 ) 


was sufficient. For the sonic boom tests, constant weights oi magnitude g- — j 25 ’ 25’ 25’ 5 5 25’ 25’ 25’ 


Collected Algorithm 

Figure C-l shows the collected deconvolution algorithm. Here, the time series is windowed using equation C-9 to 
reduce truncation smearing, Next, the transform is performed, and the spectra are deconvolved using the algorithm 
of equation C-l 7. The upper half of the spectrum is created using equation C-22, and the inverse transform is 
performed. Any Gibb's ringing in the resulting time series is eliminated using the simple filter of equation C-23. 
Finally, any effects caused by the original taper window are masked off of the deconvolved time history. 
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APPENDIX D 


SIMULATION OF TUBE DYNAMICS FOR A STEP INPUT 

To determine the linear response region of the F-16XL airdata system, it became necessary to compute the 
longitudinally and radially averaged Reynolds number in the noseboom tubing as a function ol step input size. This 
appendix describes the computational scheme used to perform the analysis. The analysis examines a small slice of 
pressure tubing where the flow is assumed to be axisymmetric. Flow perturbations are small. In addition, flow is 
parallel to the tubing walls. For axisymmetric channel flow, the Navier-Stokes equations momentum equation 
(ref. 16) may be written as 

_3[w(r, x, r)] d[p(x,f)] _ If d [~ du(r, x, Q l) n>l) 

P_ dt + ~x ^t>L r dr 1) 


where u(r, x, t) is the flow velocity in the tube, r is the radial distance from the axis of symmetry, p. is the dynamic 
viscosity, and p is the mean flow density. The longitudinal shearing stress is given by 
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Taking the radial average of equation D-3 
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where x , is the wall-shearing stress, and R is the tube radius at station j: . Defining 
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In addition, the radial average of equation D-l is written as 
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where, U{x, t) is the radial average of the velocity at station x and time t. Recalling from appendix A that the 
acoustical impedance, R a , is defined as 

_ 2 P R ey C f (A-3) 
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then 
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Across a section of tubing of length, L, if equation D-8 is integrated with respect to x, 
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Defining 
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the expression for the average velocity as a function of the pressure difference across the tubing section results 
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Recall from appendix A that for laminar conditions 
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and for turbulent conditions reference 15 gives 
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(D-l 2) and (D- 13) 


Because the average Reynolds number is a function of U (t), in the general case, equation D-l 1 is nonlinear and 
must be numerically solved to simulate the flow as a function of time, A predictor-corrector method is used to solve 
equation D-ll. The equation is first discretized at a time interval AT using explicit time differences and is used to 
predict the response based on the pressure differential and the acoustical impedance at the previous time frame 
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The Reynolds number and acoustical impedance are then re-evaluated using the predicted velocity. 
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where (ref. 1 5) 
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For 2300 < R ey < 4000, a linear interpolation of the friction coefficients was used. The result of equation D-14 is 
corrected using the new value for the acoustical impedance from equation D-16 and the pressure differential from 
the current data frame. 
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After equation D-19, the Reynolds number and acoustical impedance values were updated using the new velocity. 
Then, the cycle was repeated for the next data frame. The algorithm defined by equations D-14 and D-19 will remain 
stable as long as 
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For this analysis, the length of the tubing. L, was assumed to be length of the tubing from the input port to the digital 
pressure transducer. The effects of branches downstream of the digital transducer were ignored. Based on these lengths 
for the static and total pressure systems, the effective diameters and volumes were evaluated by iteratively solving 
equations D-14, D-15, and B-l using the experimental data from table 2 to evaluate the time constants and natural 
frequencies. The equations were solved in an ad-hoc manner until a consistent match of the system parameters resulted. 
The resulting geometry parameters are summarized in table D-l. 


Table D-l . Parameters used for simulation of tube dynamics. 


Pressure 

5 

(0 n , rad/sec 

,, . 3 

V in 

L, in. 

D eff 

Total 

0.841 

17.47 

6.8 

140 

0.114 

Static 

2.254 

8.485 

24 

140 

0.100 
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APPENDIX E 


NOMENCLATURE 


ANS astro-navigation system 

CAD central airdata 

dc direct current 

FFT fast Fourier transform 

INS inertial navigation system 

PCM pulse code modulation 


Symbols 


A[0],A[l],A[2\ 
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l LL 


l OO 
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] 
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L 
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difference equation coefficients 

tube cross-section area 

expectation operator coefficient at x = L 

expectation operator coefficient for x = L to 0 

expectation operator coefficient for x = 0 to L 

expectation operator coefficient at x = 0 

sonic velocity 

skin friction coefficient 

tube diameter 

expectation operator 

Fourier transform operator 

Gibb’s weighting constants 

optimal filter function 


cost function 
time index 
tube length 
summation index 
free stream Mach number 
data point index 
number of data points 
pressure, time domain 
pressure, frequency domain 
radial distance 
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tube radius 
acoustic impedance 
Reynolds number 
ideal gas constant 
signal-to-noise ratio 
frequency domain variable 
time 

truncation time limit 
temperature 

cosine taper truncation window 
tube velocity 

radial average of tube velocity 
volume 

longitudinal distance 

Z-transform variable 

square of damped radian frequency 

weighted averages 

ratio of specific heats for air 

overpressure, psi, psf 

sample interval 

measurement noise function, time domain 
measurement noise function, frequency domain 
tube azimuth 

truncation window degree of freedom 
weighted averages degree of freedom 
dynamic viscosity 
damping ratio 

longitudinal average of tubing density 
standard deviation 
time lag 

convolution parameter 
longitudinal shearing stress 
wall-shearing stress 
convolution function, time domain 
convolution function, frequency domain 
mean square of measured pressure, psf" 
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ft) 


frequency 
natural frequency 


Subscripts and Superscripts 

estimated parameter 

* conjugate of complex number or function 

predicted step 
longitudinal average 
eff effective 

/ parameter evaluated for flight conditions 

g parameter evaluated for ground test conditions 

i index number 

im imaginary part of complex number 

L at end of tube of length L (measured) 

r real part of complex number 

0 at beginning of tube (input) 
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